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We investigate tunneling properties of Bogoliubov phonons in a Bose-Einstein condensate. We 
find the anomalous enhancement of the quasiparticle current J q carried by Bogoliubov phonons near 
a potential barrier, due to the supply of the excess current from the condensate. This effect leads to 
the increase of quasiparticle transmission probability in the low energy region found by Kovrizhin 
et ah. We also show that the quasiparticle current twists the phase of the condensate wavefunction 
across the barrier, leading to a finite Josephson supercurrent J s through the barrier. This induced 
supercurrent flows in the opposite direction to the quasiparticle current so as to cancel out the 
enhancement of J q and conserve the total current J — J q + J s . 
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Phonon in a superfluid system is an example of Gold- 
stone mode which appears in various fields of physics 
associated with spontaneous broken symmetry It is 
a manifestation of broken (global) U(l) symmetry which 
underlies the macroscopic quantum nature of the system, 
and a key to understand the low energy properties of su- 
perfluids. In particular, in a Bose superfluid, it plays fun- 
damental roles for the superfluidity [2] . Phonon has been 
observed in various systems such as superfluid 4 He 
superconducting films [4], atomic superfluid fermi gases 
[Hj], as well as Bose-Einstein condensation (BEC) of cold 
atomic gases [ij . Due to the high degree of controllability, 
BEC of cold atomic gases offers an opportunity to study 
novel properties of phonons in the superfluid phase. 

Bogoliubov phonon [7[ in a BEC has been a long- 
standing issue of investigation in cold atomic gases [8[. 
In the last few years, quantum tunneling of Bog oliubov 
phonon has attracted much attention [9I ToL mm m 
14[ . Since Bogoliubov phonon is a collective excitation of 
a BEC, its tunneling property has specific features which 
are quite different from that of free particles. In fact, 
an anomalous tunneling property of Bogoliubov phonon 
has been predicted in [9|, LL0[ . It has been shown that the 
transmission probability of Bogoliubov phonon through 
a potential barrier increases at low energies and always 
unity in the zero-energy limit, irrespective of the height 
of the barrier 0, Hoi. Although several mechanisms 
were proposed in [l(J, LUj , the underlying physics of this 
anomalous tunneling has not been understood yet. 

In this paper, we report alternative anomalous tun- 
neling properties of Bogoliubov phonons. We find that 
the quasiparticle current is not conserved, but greatly 
enhanced near the potential barrier at low incident ener- 
gies, due to the excess current supplied from the conden- 
sate. This anomalous enhancement of the quasiparticle 
current increases the transmission probability of Bogoli- 
ubov phonon in the low energy region, which is consistent 
with the tunneling property in [9|, llOj . In addition, we 



show that the quasiparticle current twists the phase of 
the condensate wavefunction, leading to a Josephson su- 
percurrent through the barrier. The excess part of the 
quasiparticle current is canceled out by this induced su- 
percurrent, so that the total current, given by the sum of 
the quasiparticle current and supercurrent, is conserved. 

We consider a tunneling of a Bogoliubov phonon at 
T = through a potential barrier which only depends on 
x. This one-dimensional potential barrier was used in a 
recent experiment [l5| . Ignoring the motion of atoms in 
the y- and ^-direction, we can treat this model as a one- 
dimensional problem. For simplicity, we ignore effects of 
a harmonic trap. This is allowed in a box-shaped trap 
[l6[. To describe the Bose-condensed phase, we divide 
the boson field operator ip(x) into the sum of the BEC or- 
der parameter ^o(x) and the non-condensate part 5i/j(x) 
pj}, as ij)(x) = ^o(x)-\-S^(x). The condensate wavefunc- 
tion ^o(x) = obeys the static Gross-Pitaevskii 
(GP) equation, given by (setting h = 1) 
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where m, /i, and U(x) represent the mass of a boson, 
chemical potential, and a potential barrier, respectively. 
g = 47ra/m is the interaction between bosons, where a(> 
0) is the 5-wave scattering length. 

In the Bogoliubov mean-field theory, the non- 
condensate part has the form, 5^ = ^ [uj(x)&j — 

Vj(x)*Oj]. Here, dj is the creation (annihilation) 
operator of an excitation in the j-th state, satisfying the 
bosonic commutation relation [d$,dt] = Sij. Uj and Vj 
are determined by the following Bogoliubov equation: 
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where ft = -^£2 + Z7(x) + 2#|^ | 2 - r 3 is the Pauli 
matrix, and Ej is the Bogoliubov excitation spectrum. 
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We consider a simple rectangular potential barrier 
U(x) = U 0(d/2 - \x\) (U > 0). In the absence 
of supercurrent, we can safely take the condensate 
wavefunction \I/o to be real. In this case, the an- 
alytic solution of Eq. (pQ) is given in Ref. [Toj . as 
^o(x) = v^tanh -d/2)/V / 2^' + arctanh7] (|x| > 
d/2), and V (x) = ^mp/cn( ^ K 2 + f3 2 x/^, q) 
(\x\ < where cn(x,q) is the Jacobi elliptic func- 

tion, (3 = * (0)/v^, 7 = *o(d/2)/Vf^, if = 

fi/g is 
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V/? 2 + 2(/7 //i-l), and q = K/v^^2 
the condensate density far away from the barrier {x — > 
±oo), and £ = 1/ \/2mgno is the healing length. 

In the absence of the barrier, Eqs. (pQ) and ([2]) give the 
chemical potential /i = gno and excitation energy = 
yj ' e p {e p + 2gno) (where £ p = p 2 /2m). The wavefunction 
has the form v(x)) = (u p , v p )e wx 

are given by 



where u p and v p 




(3) 




FIG. 1: (Color online). Calculated transmission probability 
W (upper panel) and phase shift 5 (lower panel). We take the 
rectangular potential barrier with (d,Uo) — (£,2/i), (£,5/i), 
(£,10/i), and (4£,2/i). 



where V is the volume of the system. 
±k ±V2m 



V^ 2 + (#n ) 2 - #n describe the ordi- 
nary propagating waves in the ±x-direction. In con- 
sidering an inhomogeneous system, however, we note 
that, besides the propagating solutions in Eq. ([3]), 
Eq. ([2]) also has other two localized solutions, having 
the form of (u p: v p ) = (—6, a), where p = ±m = 

2mA 



The normalization of 
2 = -l/V. When 



these localized states is given by u 2 — v 2 
the incident Bogoliubov phonon comes from x 
the asymptotic solution is given by 
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Here, r and t are, respectively, the reflection and trans- 
mission amplitudes, satisfying |r| 2 + |t| 2 = 1. As dis- 
cussed later, this condition results from the conservation 
of energy flux. In Eq. A and B are the amplitudes 
of the localized components near the potential barrier. 

We numerically solve the Bogoliubov equation Eq. ([2]) 
so as to satisfy the asymptotic solution in Eq. (jlj). For 
this purpose, we employ the finite element method. 

Figure [1] shows the calculated transmission probability 
W = |t| 2 , as well as the phase shift S = arg(t) as func- 
tions of the incident energy E. The anomalous tunneling 
behavior discussed in [H, [HI can be clearly seen in Fig. HJ 
i.e., W -> 1 and S -> when E -> 0. Around E = 0, 
one can see the enhancement of W, the region of which 



is wider for weaker potential barrier. When the incom- 
ing energy E is very large (E ^> /i) , since the Bogoliubov 
phonon loses its collective nature, the tunneling property 
is close to that of a single particle. 

The upper panel in Fig. [1] shows the resonance tun- 
neling behavior (W = 1) at finite energies for (d, Uo) = 
(4£,2/x). At each resonance energy, u has a large ampli- 
tude in the barrier while v is suppressed. This is quite 
different from the case of anomalous tunneling at E ~ 0, 
where u and v monotonically decrease under the barrier 
approaching the zero- mode ^o- It was proposed that the 
anomalous tunneling is due to the quasiresonance scat- 
tering by the potential wells formed near the barrier • 
However, we do not find any signature of it in u and v. 

To see the low energy quasiparticle transmission in 
more detail, we directly calculate the quasiparticle cur- 
rent J q . First, we discuss current conservation for Bogoli- 
ubov phonons. In the Bogoliubov theory, the total num- 
ber density n = (V^) and current J = (1 / m)lm(^ d x ^} 
are respectively given by 



(5) 



J = J s + E ( J %- - - E • ( 6 ) 



(l/m)lm(Ujd x Uj), 
2 is the 



Here, n Uj = \uj\ 2 , n Vj = \vj\ 2 , J u 
and J Vj = (l/m)lm(Vjd x Vj), where n s = |^o 
condensate density, and J s = (l/m)Im(^ r o^^ r o) is the 
supercurrent density associated with the phase-twisted 
condensate wavefunction. The total number density n 
and total current J satisfy the continuity equation dtn + 
d x J = 0. Since the second terms in Eqs. ([5]) and (|6|) 
describe the quasiparticle contributions, the quasiparticle 



density and quasiparticle current are, respectively, given 



by n q = n u - 



and J q = J u 



J v . In a uniform system, 



the creation of a Bogoliubov phonon with momentum k 
induces the quasiparticle current J q = k/(mV). (Note 
that J u — (/c/m)a 2 , J v — (k/m)b 2 .) The last term in 
Eq. ([5]) is the so-called quantum depletion at T = 0, 
which is the non-condensate density originating from the 
repulsive interaction between bosons. 

Using the time- dependent Bogoliubov equations, 
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one obtains the continuity equations for u and v, as 
dtn u + d x J u = 5/2 and dtn v — d x J v = 5/2, where the 
source term has the form 5 = — 4glm (S&qU*v) . Thus, 
the continuity equation for quasiparticles is obtained as 



d t n q + d x J q = 5 



(8) 



In a uniform case, one finds 5 = 0, so that the number 
of quasiparticles is conserved. On the other hand, as will 
be shown later, since the source term 5 is finite near the 
barrier in our inhomogeneous problem, the number of 
quasiparticles is not conserved. In particular, in the sta- 
tionary state (dtn q = 0), Eq. (|HJ) indicates that there is a 
source supplying excess quasiparticle current in addition 
to the incoming current injected from x = — oo. 

In contrast to the non-conserved quasiparticles, the 
continuity equation for energy density p q = E(n u — n v ) 
has no source term, as dtp q + d x Q q = [10[, where 
Q q = E(J U + J v ) is the energy flux carried by quasi- 
particles. Namely, Q q is conserved in the stationary 
state. Using the asymptotic form in Eq. (|4]), we obtain 
Q q = E{k/m)(a 2 +b 2 )(l-\r\ 2 ) and J q = (k/mV)(l-\r\ 2 ) 
for x — > — oo, and Q q = E(k/m)(a 2 + b 2 )\t\ 2 and 
J q = (k/mV)\t\ 2 for x — > oo. In the stationary state, 
noting that Q q is conserved, one obtains |r| 2 + \t\ 2 = 1. 
From this result, we find that J q (x = — oo) = J q (x = oo). 

The upper panel in Fig. [2] shows the excess quasipar- 
ticle current AJ q (x) = J q (x) — J q (oo). Near the barrier, 
the quasiparticle current is enhanced, which is remark- 
able in the low energy region. This enhancement occurs 
when the source term 5 is finite, as shown in the lower 
panel in Fig. O Indeed, from Eq. dHJ, we find that AJ q is 
related to 5 as AJ q (x) = J^ OQ S{y)dy. The source term 
literally works as a source when 5 > 0, while it works as 
a drain when 5 < 0. In the barrier region (\x\ < d/2), 
the source term vanishes, so that the magnitude of the 
enhanced quasiparticle current is constant, as shown in 
the upper panel in Fig. [2l We note that the excess com- 
ponent AJ q of the quasiparticle current is completely ab- 
sorbed by the source term in the right hand side of the 
barrier. Thus, the magnitude of the outgoing quasipar- 
ticle current far away from the barrier reduces to that of 
incident quasiparticle current, as expected. In the low 
energy limit, the maximum excess quasiparticle current 




FIG. 2: Upper panel: Quasiparticle current AJ q (x) = J q (x) — 
J q (— oo) measured from the value at x = — oo. Lower panel: 
Spatial variation of the source term S. We take the rectan- 
gular barrier with (d,Uo) = (£,10//). The dash-dotted line 
indicates the potential barrier U(x). 



AJ q (a; = 0) is proportional to k. Thus, AJ q (0)/(k/mV) 
in Fig. [2] approaches a constant height which becomes 
larger for larger potential barrier. 

The transmission probability of Bogoliubov phonon 
increases when the quasiparticle current is enhanced 
near the barrier. Comparing Fig. [2] with the result for 
(d, Uo) = (£, 10/i) in Fig. [TJ one finds that the energy 
(E ~ 0.1 /i) where the high transmission probability as- 
sociated with the anomalous tunneling becomes remark- 
able coincides with the energy where the excess current 
in the barrier (normalized by k/mV) becomes large. 

The excess quasiparticle current is supplied by the 
condensate. To see this, we note that the divergence 
of Eq. (j6j) in the stationary state gives d x J = d x J s + 
Y^jj Sj(a^a.j) + \ J2j Sj. Since the supercurrent is con- 
served (d x J s = 0), it reduces to d x J = in 

the ground state ((a^aj) = 0). Since the total cur- 
rent J is conserved, the sum of the source term van- 
ishes (%2j = 0). When quasiparticles are excited 
((aj&j) 0), one obtains d x J ^ 0, which contradicts 
with the conservation of the total current. This inconsis- 
tency is actually eliminated when one includes effects of 
quasiparticles on the condensate, as 
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In this GP equation, the last term on the left side is 
the quasiparticle contribution, which originates from the 
so-called anomalous average term ^((5^^^)^q (which is 
neglected in deriving the GP equation within the Bo- 
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FIG. 3: (Color online). Phase difference (f) across the barrier 
as a function of the incident energy E. No = noV is the 
number of condensate atoms. We take Uo = 2/i, 5/i, and 10/i 
with a fixed barrier width d — £. The inset shows the phase 
6{x) for (d, Uo) = (f , 10/x) and £7 = 0.01/x. 

goliubov approximation 0, Using Eq. ©, we 

find that the continuity equation for J s is modified as 
d x Js — Sj{o^jO,j) . This result gives the expected 

conservation of the total current d x J = in the presence 
of Bogoliubov phonons. This justifies the modification of 
the GP equation shown in Eq. ([9]). 

Using Sj = d x J q j and integrating the new continuity 
equation for J s in terms of x from — oo to x, we obtain 
AJ S = — J2j AJ q j(aJ<%), where AJ s (x) is the deviation 
of supercurrent from the value based on Eq. (pp). This 
equation shows that the sum of supercurrent and quasi- 
particle current is always conserved. Thus, the excess 
quasiparticle current shown in the upper panel in Fig. 
[2] is found to be supplied from the condensate. At the 
same time, the quasiparticle current also affects the su- 
percurrent. Since the quasiparticle current is enhanced 
near the barrier, the induced supercurrent flows only near 
the barrier, in the opposite direction to the quasiparticle 
current to conserve the total current. 

If we assume that the influence of quasiparticle changes 
the condensate wavefunction by a phase shift 0(x), the 
induced supercurrent is given by AJ S = (l/m)n s (d x 0) 
[2Q( | . Thus, we find that 0(x) is given by 

6{x) = -mY J (fjy^f) <«K-> • ( 10 ) 

The inset in Fig. [3] shows the spatial variation of Nq9(x), 
where N = n^V is the number of condensate atoms. 
Here, we assume that only one Bogoliubov phonon with 
energy E exists. The phase 9{x) sharply varies around 
the barrier. This clearly shows that the Bogoliubov 
phonon twists the relative phase of the condensates on 
the left and right sides of the barrier. Since their cou- 
pling is weak at the barrier region, the relative phase is 
subject to change by the tunneling quasiparticle. The in- 
duced supercurrent A J s can be naturally regarded as the 
Josephson current originating from the phase difference 
(j) = 0(x = —00) — 0(x = 00). 



Figure [3] shows the normalized phase difference 
7V as a function of E. It is enhanced at low 

energies similarly to the normalized excess quasiparticle 
current AJ q /(k/mV) in Fig. [2] and indeed, normalized 
supercurrent AJ s /(k/mV). Actually, the maximum su- 
percurrent AJ s (x = 0) is proportional to <p at low ener- 
gies. This result justifies regarding J s as the Josephson 
current due to (j). In addition, Fig. [3] shows that, at the 
same energy, (j) is larger for higher barriers. This implies 
that, since the coupling of the condensates across the bar- 
rier becomes weaker as the barrier gets higher, the phase 
of the condensate wavefunction is easily twisted by the 
Bogoliubov phonon for higher barriers. 

In summary, we have investigated tunneling effect of 
Bogoliubov phonon in a BEC at T = 0. We found anoma- 
lous tunneling properties of Bogoliubov phonons. The 
quasiparticle current of Bogoliubov phonon is enhanced 
near the potential barrier due to the supply from the 
condensate. This enhancement is remarkable at low inci- 
dent energies, and explains the increase of the transmis- 
sion probability of Bogoliubov phonon in [9|, [l0[. Fur- 
thermore, we have shown that the quasiparticle current 
twists the phase of the condensate wavefunction across 
the barrier, which induces the counterflow of Josephson 
supercurrent to conserve the total current. This twisted 
phase could be observed in the interference pattern when 
a BEC is released from the trap potential as in [211 ]. 
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